Global controllability with a single local actuator 



S. G. Schirmer/'j^] I. C. H. PuUen,^ and P. J. Pemberton-Ross^ 

' Dept of Applied Maths and Theoretical Physics, 
University of Cambridge, Wilberforce Rd, Cambridge, CBS OWA, UK. 
(Dated: June 6, 2008) 

We show that we can achieve global density-operator controllability for most TV-dimensional bilin- 
ear Hamiltonian control systems with general fixed couplings using a single, locally- acting actuator 
that modulates one energy-level transition. Controllability depends upon the position of the actua- 
tor and relies on the absence of either decompositions into non-interacting subgroups or symmetries 
restricting the dynamics to a subgroup of 6il(A''). These results are applied to spin-chain systems 
and used to explicitly construct control sequences for a single binary-valued switch actuator. 
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I. INTRODUCTION 

The ability to control the dynamics of quantum sys- 
tems is a long established objective in areas as diverse as 
molecular chemistry and quantum computing among oth- 
ers. Control in practice comprises various related tasks 
such as transforming a system from a given initial state 
to a desired target state, implementing a desired uni- 
tary operator, or optimising the expectation value of a 
selected observable. The manner in which control is ef- 
fected depends on the system but a common approach 
for quantum systems is the application of external elec- 
tromagnetic fields. In the diabatic control regime these 
drive transitions between different states of the system, 
and control can be achieved by adjusting the amplitude 
and phase of the field(s) as a function of time in a way 
that maximizes constructive interference of various ex- 
citation pathways that lead to a desired outcome, while 
maximizing destructive interference for all others. 

Although the ultimate goal of control is usually to find 
a control field that steers the system in the manner re- 
quired to achieve the objective, the question of what tasks 
can be accomplished for a given system with a given set 
of actuators, is of fundamental interest. A key concept in 
this regard is that of controllability. A substantial num- 
ber of papers have been devoted to studying this issue for 
both classical and quantum systems, establishing various 
notions of controllability and general algebraic criteria for 
them, and showing that particular types of systems are 
controllable [1] El E HE [i E [8] . On the latter front, it 
has been shown, for instance, that any system with dis- 
tinct transition frequencies and a connected transition 
graph is controllable [SJ [TU] . For an n- level system this 
requires at least n—1 transitions with non-zero probabil- 
ities. It has also been shown that these requirements can 
be relaxed in many cases jllj , and more recently indirect 
controllability has been studied IT2]. 

One remaining area of interest is global controllability 
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with a small number of local actuators. A motivation for 
this type of scenario could be a chain or array of quan- 
tum dots with control electrodes to locally manipulate 
the dynamics of one or a few quantum dots. In the ideal 
case, one might consider separate control electrodes for 
each quantum dot, as well as separate electrodes to mod- 
ulate all the interactions between pairs of adjacent quan- 
tum dots, as proposed by Kane in |13j and many other 
quantum computing architectures since. Leaving aside 
the often considerable challenge of finding optimal con- 
trol schemes and fighting decoherence, with sufficiently 
many local actuators almost any (Hamiltonian) quantum 
system is controllable, at least in principle. However, in 
many cases it is impractical or even impossible to have 
a large number of individual local actuators such as con- 
trol electrodes. Rather, one would like to make do with as 
few local actuators as possible to simplify the engineering 
design and reduce deleterious efi^ects such as decoherence 
and crosstalk, for example. 

Motivated by this problem we investigate the question 
of controllability of a finite-dimensional model system 
with the smallest number of simple actuators whose effect 
is strictly confined a local perturbation of the Hamilto- 
nian. We also note here that by local we mean localized 
in space, affecting a single transition, for instance, not si- 
multaneous local operations on many individual elements 
such as qubits as is common in global control schemes. 
We effectively show that in most cases a single local ac- 
tuator is sufficient to ensure controllability of the system 
as a whole, provided the latter is not decomposable into 
non-interacting parts, and does not exhibit dynamical 
symmetries that its evolution to a subgroup of the uni- 
tary group (IX(Af) or &'d{N)). Many systems with fixed 
interactions connecting its parts such as chains of quan- 
tum dots etc with fixed non-zero couplings between ad- 
jacent dots satisfy this connectedness requirement, and 
the disorder present in most realistic systems is likely 
to ensure that there are no special dynamical symme- 
tries to worry about in most cases. For these systems 
our controllability analysis suggests that the entire sys- 
tem can be controlled by modulating a single transition 
with a local actuator. Although the explicit controlla- 
bility proofs given apply to specific model systems, the 
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same arguments are applicable to many other model sys- 
tems, suggesting that a large class of systems with fixed 
couplings may be controllable using a very small number 
of fixed local actuators. We conclude with an explicit ex- 
ample of constructive control with a single binary switch 
actuator. 



II. DEFINITIONS AND BASIC RESULTS ON 
CONTROLLABILITY 

We restrict ourselves here to control problems that 
can be classified as open-loop Hamiltonian engineering 
problems and systems subject to Hamiltonian dynam- 
ics. Open-loop control engineering means that we aim to 
design control fields relying only on (presumed) knowl- 
edge of the initial state of the system and the dynamic 
laws governing its evolution in the presence of the con- 
trol fields, without any feedback from measurements. We 
furthermore assume the state space of the system is a 
finite-dimensional Hilbert space, H ~ C^. The state of 
the system in this case can be represented by a density 
operator p, i.e., a positive unit-trace operator acting on 
TL, and its evolution is governed by the quantum Liouville 
equation 



(1) 



where [A, B] — AB — BA is the usual matrix commuta- 
tor and Cd — Q for a Hamiltonian control system. The 
operator H\f{t)\ is the total Hamiltonian of the system 
subject to the control fields f(t). For control-linear sys- 
tems we have the perturbative expansion 



H[m] = ^0 
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where Ha is the internal Hamiltonian of the system and 
Hm, m > 0, are the interaction terms. 

Hamiltonian dynamics constrains the evolution of den- 
sity operators p{t) to isospectral fiows 



p{t) = U{t,t„)poU{t,to) 



(3) 



since the evolution operator U{t,to) must satisfy the re- 
lated Schrodinger equation 



ih--Uit,to) ^ H[{{t)]Uit,to) 
at 
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and is hence restricted to the unitary group !d{N) . Due to 
this fundamental restriction it is clear that the maximum 
degree of state control we can achieve for this system 
is the ability to interconvert density operators with the 
same spectrum, which is achieved if we can implement 
any unitary operator in the special unitary group &ii{N) 
of unitary operators with determinant 1 as abelian factors 
do not affect the isospectral flow. It is also not diflicult to 



show that any proper subgroup of &il{N) is not sufficient 
to interconvert any two generic density operators with 
the same spectrum. 

To properly define the notion of controllability we need 
some concepts from Lie group / algebra theory. A Lie 
algebra is a vector space over a field endowed with a bi- 
linear composition [x, y] that satisfies the Jacobi identity 



'[x,y],z] + [[y,z],z] + [[z,x],y] = 0. 



It is easy to see that the anti-Hermitian matrices iHo and 
iHi generate a Lie algebra £ which must be a subalge- 
bra of the Lie algebra of skew-hermitian matrices u(iV), 
and if iHo ^md iHi have zero trace, £ will be a subal- 
gebra of the trace-zero, anti-Hermitian matrices su(iV), 
which can be regarded as the tangent space to the Lie 
group 6il{N) at the identity via the exponential map 
X G su{N) 1-^ exp(a;) G 6ii{N). Therefore, we can argue 
that if the iHg and iHi — or their trace-zero counterparts 
Hm = Hm — Ti{Hm)lN — generate the entire Lie al- 
gebra su{N) then we can in principle dynamically gen- 
erate any matrix U E &ii{N). Hence, a system is said 
to be density matrix controllable or simply controllable 
if the Lie algebra generated by iHo and iHi is su(iV). 
These Lie algebraic criteria are useful as they are easy to 
check by quite straightforward calculations. In principle, 
these can be done numerically for a given set of Hamilto- 
nians but for higher dimensional systems the calculations 
can be time-consuming and the accuracy very limited. It 
is therefore desirable to have more explicit criteria that 
guarantee controllability for certain classes of systems, 
and several such results exist. 

For example, consider a simple finite-dimensional sys- 
tem (dimTi = N) with a control-linear Hamiltonian of 
the form Hq + f{t)Hi. Choose a basis such that Ho is 
diagonal with energy levels _E„, n = 1, . . . , N, and transi- 
tion frequencies LUmn — En—Em- If ^^o is regular, i.e., has 
non-degenerate eigenvalues then we can associate each 
1-dimensional eigenspace with the vertex of a graph and 
interpret the non-zero elements in the matrix represen- 
tation of the interaction Hamiltonian Hi (with respect 
to the eigenbasis of Hq) as edges of a transition graph. 
In this case a sufficient condition for controllability is 
that Hq be strongly regular, i.e., have distinct transition 
frequencies w„i„ 7^ tOm'n' unless {m,n) = {m',n'), and 
the transition graph as defined above be connected [9j. 
The conditions of uniqueness of the transition frequencies 
can be slightly relaxed in that we only need to consider 
the transition frequencies of those transitions that occur 
with non-zero probability. This is a useful result as it is 
very easy to check, although it is important to remember 
that it provides only a sufficient, not a necessary condi- 
tion, and indeed many systems that do not satisfy these 
conditions are controllable. For example, given a system 
with a Hamiltonian of the form H[f{t)\ = Ho + f{t)Hi, 



where 



N 



^0 - ^^^n|«)(?^|, (5a) 

ri=l 

^fi - I] c^n[|«>("+l| + l"+l)(«|], (5b) 



(a) Global control field simultaneously driving all transitions 



the graph connectivity result allows us to conclude that 
the system is controllable provided the energy levels of 
the system are such that the frequencies of all transi- 
tions between adjacent states are distinct and dn ^ for 
n=l,2,...,A^— 1. In principle this controllability result 
can be explained in terms of frequency-selective control. 
If all the possible transitions have different frequencies 
then we can imagine a field resonant with a particular 
transition frequency as selectively driving only the reso- 
nant transition. Thus, we can implement all the genera- 
tors Xn.n+i and yn,n+i of the Lie algebra and therefore 
the entire Lie algebra, implying controllability. However, 
it has been shown that many systems that do not satisfy 
this condition such as the truncated harmonic oscillator 
with dn = \[n are controllable despite all allowed tran- 
sitions having the same frequency [TT]. What permits 
controllability in this case, despite the lack of any fre- 
quency selectivity, are the differences in the transition 
strengths d„- If both the frequencies and strengths are 
the same for all transitions, or even if they satisfy certain 
not necessarily obvious symmetries, then controllability 
is indeed lost [14 . 

Furthermore, both the graph connectivity and the ex- 
plicit Lie algebraic results mentioned above apply only to 
systems where the field drives all the possible transitions, 
i.e., is global as shown in Fig.[Ta| Given a particular con- 
trol field that drives only a single transition, or a subset 
of all possible transitions (see Fig. lb I, we cannot draw 



any conclusions even if the transition graph of the system 
appears connected and the transition frequencies are dis- 
tinct. Of course, if we have many local control fields, each 
selectively driving a single transition, as shown in Fig Ic 
then it is again obvious that the system is controllable. 
It is not obvious, however, under which conditions the 
ability to control a single transition of a large connected 



system as in Fig lb is sufficient for controUabilty. 



III. 



CONTROLLABILITY FOR SINGLE LOCAL 
ACTUATOR 



In this section we consider a model system with N dis- 
tinct states that are permanently coupled in some form, 
such as an array of quantum dots. For our model we first 
assume coupling of nearest neighbour type, which leads 
to a drift Hamiltonian of the form 



(6) 



(b) Single local actuator driving a single transition 



(c) Many local actuators driving individual transitions 

FIG. 1: Schematic of A^-state system with coupling 
between adjacent states with (a) global control field, (b) 
a single local actuator and (c) many local actuators. 



with i7o and ii\ as in Eq. ([s]) , and a single local actuator 
modulating the coupling between states |r) and |r-|- 1), 



Ar = |r)(r+l| + |r + l)(r|. 



(7) 



Thus, for a single local actuator positioned between r 
and r + 1 we have the total Hamiltonian 



(8) 



Example 1. The Hamiltonian of the first excitation sub- 
space of a spin chain of length N with nearest neighbour 
coupling of isotropic Heisenberg form given by the cou- 
pling constants d„ > for n = I, . . . , N — 1 is 



N 



N-1 



^0 = ^Sn|")(«|+ ^ dn[\n){n+l\ + \n+l){n\], (9) 



71=1 



where the energy levels are explicitly 

1 NT^ , 1 



E„ 



di--{d„ 



dn) 



(10) 



and we set do — d^ ~ 0. Assuming we have a local ac- 
tuator that allows us to modulate the coupling between 
spins r and r -|- 1, the total Hamiltonian is of the form 
Eq. ^ with Ar of the form Eq. The first excitation 
subspace Hamiltonian for spin chains with dipole-dipole 
interactions is also of form ^ but with different energy 
levels. Similar results hold for any spin chain decompos- 
able into excitation subspaces. 

Again, it is quite obvious that TV — 1 independent local 
actuators of this type, controlling the coupling between 
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spins n and n+1 in the chain, will suffice for the system to 
be controllable, but in fact, a single such actuator suffices 
in most cases. 

Theorem 1. A quantum system with Hamiltonian 
^[/(^)] = ^0 + f{t)^r with Aq and Ar as above is con- 
trollable if uJr 0, dn and rf^+i 7^ c^r-i- 

Proof. We show that the trace-zero anti-Hermitian ma- 
trices iVo and iVr defined by 



Tr(Ao) TvjA^) 

1^0—^0 ^N, Vi — Ar — 1. 



N 



generate the Lie algebra £ — su{N). To this end it suf- 
fices to show that the Lie algebra £ contains the 2{n — 1) 
generators a;„ = Xn,n+i and ?/„ = y„,„+i of su{N), where 
the basis elements of su(A^) are defined as usual, 

Xmn = \n){m\ ~ \m)(n\, 
Vrnn = i{\n) {■m\ + \m){n\), 
hn ^ \n){n\ - |n+ 1|, 

for 1 < m < n < iV. Let V^°^ = i{Vo - drVi). We have 
iVi — Ur ^ ^ and 

-'^O = [Ur, Vq^''] = dr-lXr-l^r+l — dr+lXr,r+2 — l^rXr 
Yq = [-'^O, Vr] = dr-lVr-l + d^+lj/r+l - 2^^^^ 

X'q = [Yo^Ur] = —dr-lXr-l^r+1 + dr+lXr^r+2 + 4,0JrXr 
Yq = [-'^Oi Vr] = -dr-lUr-l - dr+ll/r+l + BuJrK 

yields Xr = {?iUJr)^'^ {Xt) + X'^) G £ and hr — 2^^[xr, Vr] G 
£ as 7^ 0. Next setting 

Yl = 3"^(4ro + Yq) = dr-lUr-l + dr+lUr+l 
Xi = [[Xr, Yi],l/r] = rfr-lXr-l 

Zi = 2-i[^i,yi] - dl_^K-x + d'+i/ir+i 
Y[ = 2-'[Z,,X,] = dl_^Vr-l + dl+^Vr+1 
X[ = 2-\Yi, Zi] = dl_^Xr-l + dl+^Xr+l, 

and ci = rf^-i — '^r+i leads to 

Yi - rfr+l^l = dr-lCiVr-l, 





-dl^ 


iXi 


= dr-lCiXr-l, 


Yi 


-dl 


-lYi 


= -dr+lCiyr+1 


x[ 


-dl_ 


iXi 


— —dj,+ lCiXr+l 



Since dj.±i 7^ 0, ci 7^ by hypothesis, we have yr±i, 
Xr±i, and hr±i — 2^^[xr±i,yr±i] in £■ Next note that 



E 



dnUn 



where is the index set {1, . . . , A^— 1} minus the subset 
{r — 1, r, r -|- 1} and we have 



Y2 = [[Zi,V^^''],Zi] = dr-2dt-iyr-2 + dr+idr+2yr+2 
V^^^ =V^'^ -d;\Y^ = tHo+ J2 dnyn + Cr+2Vr+2, 



with /(^) the index set /^""^^ minus {r — 2,r -I- 2} and 

Cr+2 - - '-^^ 



c?r+2(l — d^+i/^r-i)- Hence 



X. 



dl 



■r+l'^r+2^r+2: 



Y2 = [X2,Zi] = dr+iCr+2Vr+2 

shows Xr+2-, yr+2 and hr+2 = 2^^[x^^'^ ,yr+2] in £. Set- 

(3) (2) 

ting Vq — Vq dr+2yr+2 HOW shows that 

yr+3 = hr+2] e £, 

hr+z = 2"^[a;r+3,2/r+3] e £. 

Repeating this procedure with Vq^^'^'^ = V^^^ — dr+kVr+k 
we obtain 



d-l^^,[K+u,vt+'^]€Z 

[Xr+k+li hr+k] G £ 
2"Ma^r+fe+l,yr+A;+l] G £ 



yr+fc+l 
hr+k+l 

for 3<fc<A^ — r — 2. To show that the elements Xr-k^ 
yr-k for 2 < k < r — 1 are in £, we note that 

yr-2 = rfr-2^7-l(^2 - dl^ldr+2yr+2) G £ 
a;,._2 = [/lr-l,yr-2] G £ 
h.r-2 2"^[a;r-2,yr-2] e £ 



and setting W^^^ = 



and VFr 



(fe+i) 







= W, 



(k) 



dr-k-iyr-k-l shows 

yr-k-1 — [Xr-k-l,hr-k] G £ 
/lr-fe-1 = 2 ^[Xr-fc-l, ?/r-fc-l] G £ 

for 2 < fc < r — 2, as desired. 



□ 



For a Heisenberg spin chain dn > for n = 1, . . . , iV — 1 
and Eq. (10 1 show that lj„ = dn-i — dn+i- Thus ujr =/= 



is equivalent to c?r+i 7^ c?r -i and we have the following 

Corollary 1. The first excitation subspace of a Heisen- 
berg spin chain of length N with coupling constants dn 
is controllable with single local actuator between spins r 
and r + I if d^+i 7^ dr^i- 

A Heisenberg spin chain with non-uniform couplings 
almost certainly satisfies c?^_x ^ c'r+i fo'" any r between 
1 and iV — 1. A chain with uniform coupling d„ = d, 
n — 1, . . . , N — 1, do = d^ = 0, satisfies this condition 
only if the actuator is placed near the end of the chain, 
i.e., r = 1 or r = — 1. However, we can generalize the 
previous theorem. 

Theorem 2. A quantum system with Hamiltonian 
^[/(^)] — ^0 + fit)Ar with Aq and A^ as above is con- 
trollable if UJr 7^ 0, dn ^ and rf^-fc-i ^r+fe+i some 
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Proof. For k = 0, i.e., if d^-i (^r+ii the result follows 
from Thmjl] If dl_i — d^^j^ we begin as in the proof of 
Thni[l]to conclude that yr G Xr — (3a;,.)^^(Xo+Xg) e 
£ and hp = 2^^[xr, Vr] G and set 



(0) _ 



^(0) 
(1) 



(0 



(0) 



dnVn 



where I^^^ is the index set {1, . . . , A^— 1} minus the subset 
{r — 1, r, r + 1}. 

Setting d'^j^j = '^r-j ior j — 1 . . . k — 1 and observing 
that we cannot separate the r + 1 to r + fc and r — 1 to 
r — k terms, respectively, at this stage we continue along 
similar lines by iterating the following set of recurrence 
relations for j = 1, . . . ,k — 1 



zf^^d-l^z^^^ 

~ hr-j + hr+j as d^.j^j = dl_j ^ 

— d/p— J dip— J ~ ']^x — J — \ .T — J -)~ 1 dip — -jCiJ'p—jh'p — J 
— dj^^jdj^-^j-^iX-p-^j ,p-\-j-\-2 ^ d^-^jUJ'p-^jh^^j 



— dj._jd.p—j—iyr—j—i 2d^_jUj'p_jhj^_j 
d.^_^j df 

(2) - ^-2 vd) 



^r+jdr+j + lUr+j + l ~ '^d^+j^r+jh-r+j 



— dfp — J ]_ — J — X " ^^ij-^'p — J I^'p — J 

+ dr+j + lUr+j + l — 2uJr+jhr+j 

xfl,^\zf\Y^] 

— dr^j-lXr-j^l + dr+j + lXr+j + l 
^(0) _ rT(^(0) 7(1)1 

— dr-j-lUr-j-l + dr+j + l2/r+j + l 
7(0) _ r,-lry(0) T^(0)i 

T/(J + 1) _ T/(J) _ 

''O — ''^0 + 1 

= -iHo - ^ dnVn 



where is the index set /'^■'^ with the subset {r — j 



1, r + j + 1} removed. Since d^^_f,_^ ^ ^l+k+i ^'^'^ 



Y, 



(0) 

fe 

(0) 
.(0) 



d^ — kXr—k ^" dj'-]-kXf-^k 
dr-kUr-k + dr+kVr+k 



(fe) _ 



where I^*'-' is the index set {1, . . . , r — fc, r + fc, . . . iV}. 

To complete the proof by showing that yr±(fe+i)i 
Xr±(k+i) and /ir-i-(j.^i') are in £, we calculate the com- 
mutators 

xi''^[Yr,Vn 

dq^_'j^dj'—]^—\Xj^_]^_\_q^_'j^j^\ d>p_f^uj^^f^h^^f^ 

— dr^kdr-\-k-\-lXr-\-k,r-\-k-\-2 ~ dr^k^r+khr^k 

y!^' ^[x^\yI^^] 

— dl_i^dr-k-iyr-k-l — 2dl_f.UJr-khr-k 
+ d^^i^dr-^k+lUr+k+l ^ '^d^j^^.LOr+khr+k 

4^' - d^ikZ^' 

= hr—k + hr+k as d^—k ~ ^r+k 



4''^[zi'\Y^^^] 

^ d-r-k'^r-k-l^r-k-l + d^^i^dr+k+l^r+k+l 
= dr-k-lXr-k-l + dr+k+l^r+k+l 

y!P^\x^,z^^] 

= dr-k-lVr-k-l + ^r+fe+lZ/r+fe+l 
= df._k_ihr-k-l + d^r+k+l^r+k + l 
= dl-k-lVr-k-l + rf?+fe+l2/r+fe+l 



which gives 

yr±(fc+l) 



(3) 



- "r-T(fc+l)^fe 



(2) 



'^r±(fc+l)(rfr±(fe+l) ^ ^rT(fe+l)) 
Xr±(k+1) = br±(fe+l), ^^■']/(2rfr±(fc+l)) 
^r±(fe+l) = [2^r±(fe+l), 2/r±(fc+l)]/2 

showing that these generators are in £. To show that the 
generators Xr±j, yr±j, and hr±j are in £ for j = 1, . . . , /c 
we set 

,^(fe+i) _ -.^(fc) , , 

V'O — ~ "r-fc-lj/r-fc-l — Ur+k+lVr+k+l 
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where /C^+i) is the index set I^''^ with the subset {r — 
k — l,r + k + 1} removed, and note that 

Xr±j = rf7±j[[yr±(j + l),^j"^],2/r±0 + l)] 
yr±] = [Xr±(j+1), [Xr±j,yr±{j+1)]] 

Finally, we show that the generators Xr+k+j, Ur+k+j and 
hr+k+j are in £ for j = 2 . . . iV — r — fc — 1 , by iterating 
the following set of recurrence relations for j = 2 . . . iV — 
r — k — 1: 

Xr+k+j = dr+k+'j[^r+k+j-l-,V^^^^ 
Vr+k+j = [Xr+k+j, hr+k+j-l] 
hr+k+j = 2 ^[Xr+k+j, Vr+k+j] 
^0 ^ ^0 ^ Ur+k+jUr+k+j 

Similarly, we show that the elements Xr-k-j, Ur-k-j and 
hr^k~j are in C for j = 2,...,r — fc — 1, by setting 
Wq'^ — Vq'^'^^^ and iterating the following recurrence 
relations for j=:2...r — fc— 1: 

Xr-k-j d;\_^[hr-k-j + l,W^''^^ + ^'>] 
Vr—k—j ^ ['^r — k—j-)^r—k—j-{-l\ 

W^'-''> = W^'-' + '''~dr^k-0r-k^, 

We have now shown that Xj , yj E £, for 7 = 1 ... — 1 
as desired, completing the proof. □ 

We note that dl_^_^ 7^ rf^+zj+i for some integer k is 
always satisfied if the system dimension is odd, N = 
2£ + 1, no matter where we place the actuator. If = 
2£ then d'^_j._i = c?r+fe+i ^o^' ^ possible only if 
r = i, i.e., if the actuator is placed in the middle, and 
the coupling constants are symmetric around the centre, 
rf|_j. = for all k. For a spin chain with strictly 

isotropic Heisenberg interaction the requirement Wj, 7^ 
is still a problem if dr-i = dr+i but controllability could 
be restored by engineering a local perturbation of the 
energy levels in the vicinity of the actuator, which may 
indeed achieved by the actuator itself. 

We can interpret these results in terms of transition 
graphs. Given a system with a tridiagonal drift Hamil- 
tonian Hq with respect to some Hilbert space basis 
{\n) : 1, . . . , N}, we can define a transition graph as be- 
fore by taking the N basis states as vertices and adding 
edges for each non-zero transition. Since the Hamilto- 
nian is tridiagonal the resulting graph is either a linear 
chain or disconnected. Connectedness is a necessary con- 
dition, and the results above guarantee controllability in 
the following cases. 

• If the chain is connected and has odd length then 
the actuator can be placed anywhere provided the 
vertices associated with the controlled edge have 
different energy levels. 




FIG. 2: Controllability of a chain with single actuator: 
model systems (a) and (d) are controllable, (b) is 
controllable if ^ E4, while (c) is not controllable 
even if E3 ^ E4 if chain has reflection symmetry, i.e., 
\d2i-k \ = \dk\ and E2i-k+i = Ek for all k as the 
associate dynamical Lie group has symplectic symmetry. 



• If the chain is connected and has even length then 
we must ensure in addition that the system does not 
admit symplectic symmetry, which is guaranteed as 
long as the actuator is not placed precisely in the 
middle of the chain. 

• For a uniform chain for which the energy levels in 
the interior of the chain are always degenerate, con- 
trollability is ensured by placing the actuator near 
cither end of the chain. 

It is worth pointing out here that the transition frequen- 
cies of the system need not be distinct. In fact, the most 
of the energy levels can be degenerate as is usually the 
case for uniform spin chains. We only require that the 
vertices of the controlled transition have different energy 
levels, which can generally be achieved by placing the 
actuator near the end of chain. 



IV. CONSTRUCTIVE CONTROL WITH 
SINGLE BINARY SWITCH ACTUATOR 

The results of the previous section suggest that a sin- 
gle actuator is often sufficient to achive the same degree 
of controllability that is achievable with many local ac- 
tuators. This result is not too surprising on purely Lie 
algebraic grounds considering that two randomly choosen 
Hermitian N x N matrices, generically, will generate the 
entire Lie algebra u(A^). Of course, the Hamiltonian ma- 
trices in our model system are far from random, and our 
Lie algebra calculations indeed show that for certain sys- 
tems such as a spin chain with uniform, isotropic Heisen- 
berg coupling between adjacent spins, controllability de- 
pends on the type of interaction and the placement of 
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the actuator. Nonetheless, that a single local actuator in 
many cases results in the same degree of controllability 
than, say, 2N — 1 local actuators to individually control 
all of the energy levels and transitions, is rather surpris- 
ing when one considers the substantially reduced control 
that such an actuator affords, and it begs the question 
whether it is possible to control such a system construc- 
tively, i.e., whether we can find a local control field /(t) 
that achieves the desired global system dynamics, and 
if there exists a solution, whether it can be practically 
implemented. 

The type of control functions that are feasible gener- 
ally depends on the specifics of the system and actua- 
tor. For laser-controlled quantum dots, for example, the 
availability of pulse shaping technology and the demon- 
strated superiority of shaped pulses over simple pulses 
in certain settings, suggests optimization routines de- 
signed to find an optimal time-dependent pulse shape 
/(t), and many such algorithms based on gradients and 
variational techniques have been proposed (see for ex- 
ample [iSl Uni [17] ) . For many other systems, especially 
voltage gate controlled systems, however, it is generally 
difficult to implement complicated time-varying poten- 
tials, and simple, piecewise constant controls that can 
be approximated by square pulses are preferable. In the 
following we consider the simplest type of such an actu- 
ator, a binary switch that switches the voltage between 
two possible values, corresponding to two fixed Hamilto- 
nians 



(Ua) 
(lib) 



Given a sequence of switching times t = (ti, . . . ^Ik) the 
corresponding evolution of the system is given by 

U{t) = U^^\h)U^^\t2) . . . U^^\tK-l)U^^\tK) (12) 

where C/("H*fe) = exp{~itkH^"^^) for to = 1,2. The 
control task in this case is reduced to find the switching 
times t to accomplish a desired task. Although analytical 
expressions for the optimal switching times are generally 
very difficult to obtain for all but very simple systems, 
numerical optimization techniques can be used to find 
suitable controls, and we have found them to be surpris- 
ingly effective in many cases. 

As a specific example, we consider the first excitation 
subspace of a spin chain of length four with a single bi- 
nary switch actuator placed between spins one and two, 
i.e., r = 1. This system is controllable with a single ac- 
tuator at r = 1 according to Theorem [T] To show that 
we can constructively control this system with a single 
binary switch actuator, we find switching time sequences 
for a complete set of generators of 6il(4). Interpreting 
the first excitation subspace of the chain as a two-qubit 
system by setting 

|0) = |00), |1) = |01), |2) = |10), |3) = |11), 



we show that it is possible to find vectors t^'*^ such that 
||C/^"^ -C/(t(^))|| < 10-4 (13) 
for the following set of six target operators 



U!j."^ e {1(E) I, Had(^/, T(g)I, /®Had, I(g)T, CNOT}, 

(14) 

where / is identity operator on a single two-level sub- 
space (qubit), T = exp(— ztt/Sctz) is a n/8 phase gate, 
and Had and CNOT are the Hadamard and CNOT gate, 
respectively, 

Ha..-L(; -1), CNOT.e-/.(;^^0). 

with and az being the usual Pauli matrices 



1 
10/' 



1 

-1 



The set of target operators ( 14 1 was chosen because it is a 



universal set of elementary gates in that any other 6il(4) 
operator can be constructed from these elementary gates, 
and the ability to implement a universal set of gates for 
SU(4) is equivalent to (density operator) controllability 
of the system. 

Table |l] shows the time vectors t*^*), as well as the 
gate operation times T = X^fcLi ^^'^ S^^^ errors as 
defined above for a system with uniform isotropic Heisen- 
berg coupling. The optimization was performed using a 
Nelder-Mead downhill simplex algorithm |TH] with mul- 
tiple initial simplices. The table shows that it is possi- 
ble to implement all of the six elementary gates with a 
fidelity > 99.99% with no more than 20 switches of a sin- 
gle on-off switch actuator in approximately 40 time units 
each, a surprisingly good result considering the minimal 
nature of the available control. Since solutions are obvi- 
ously not unique, and the minimum control time or num- 
ber of switches required to achieve the control objective 
are unknown, even better solutions probably exist. The 
non-uniqueness of the solutions can be exploited to sat- 
isfy additional constraints such as minimum pulse lengths 
(switching cannot be arbitrarily fast) etc. 



V. CONCLUSION 

We have shown that a certain class of systems of 
Hilbert space dimension N is controllable with a single 
local actuator. In particular, the results show that it 
is usually not necessary to be able to control all transi- 
tions, and a single local actuator in fact suffices in most 
cases to achieve controllability. That is to say, we do not 
require N — 1 or more independent local actuators, or 
global actuators acting on the entire system. 

The results establish theoretical minimum require- 
ments for controllability for a class of systems that in- 
cludes many types of spin chains and other systems with 
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I®/ 


Had (g) / 


T® J 


J® Had 


/®r 


CNOT 


error 


6.02407e-05 


9.93462e-06 


9.41944e-10 8.86637e-07 


1.10773e-06 


1.31773e-06 


duration 


40.5351 


37.9537 


41.166 


41.1328 


42.5368 


39.3569 


ti 


0.731996 


3.94518 


1.79446 


3.08601 


1.30764 


3.34576 


t2 


2.03884 


2.20021 


1.79932 


3.13305 


1.1069 


0.0179813 


ts 


3.52727 


0.0384191 


0.0730935 


0.701478 


0.518925 


2.59171 


U 


1.38628 


1.07432C-07 


1.71885 


3.62498 


5.85085 


3.23448 


t5 


3.39919 


0.680856 


2.07051 


2.45712 


0.396729 


1.46693 


te 


0.951534 


3.04816 


0.747468 


0.68558 


7.37392 


0.212212 




1.35113 


1.292 


1.84047 


1.3746 


1.13031 


5.0851 


ta 


0.575672 


1.86256 


2.53341 


1.12801 


1.16712 


3.07975 


t9 


3.38307 


4.14879 


4.73792 


3.50997 


0.802765 


2.75667 


tio 


0.0365974 


0.356856 


1.3432 


1.92944 


4.08279 


0.439889 


tu 


3.62131 


1.02202 


1.39084 


5.57909 


1.27132 


3.25423 


tl2 


0.93505 


0.0453206 


0.320722 


0.298252 


2.70023 


2.41685 


tl3 


1.75377 


2.13701 


4.15595 


0.987279 


4.67647 


1.04768 


tl4 


5.19515 


1.24291 


0.533115 


0.26934 


0.705919 


1.31426 


tl5 


4.5099 


0.101593 


1.03574 


1.7998 


1.01477 


2.6859 


tie 


1.01899 


4.40131 


7.58673 


4.66334 


1.78438 


0.732592 


tir 


4.01314 


1.07241 


4.77061 


0.135612 


1.02283 


0.16703 


tl8 


0.991019 


5.83516 


0.857316 


1.31499 


1.02426 


0.770284 


tig 


0.705316 


1.6229 


1.7735 


3.38444 


2.16687 


2.32287 


t20 


0.409887 


2.89999 


0.082803 


1.07044 


2.43177 


2.41471 



TABLE I: Gate errors (1— gate fidelity), total time T required to implement respective gates, and vector of switching 
times tk to implement a universal set of elementary gates with 20 switches for (the first excitation subspace of) a 

uniform isotropic Heisenberg spin chain of length four. 



non-trivial fixed couplings between adjacent elements. 
As systems with fixed interactions are generally much 
easier to engineer than systems with individually tunable 
transitions, this is a promising result. 

Although the controllability proof is an existence 
proof, we have further demonstrated that is is possible 
to constructively control a system with a local actuator, 
even if the actuator is limited to a binary switch, for a 
four-level system, where we have shown that it is possi- 
ble to implement a complete set of generators of 611(4) 
with fidelities > 99.99% using a single binary switch ac- 
tuator, with no more than 20 switches per gate required. 
Although it would be desirable to have analytic expres- 
sions for the switching times, it appears that numerical 
optimization techniques are quite effective in finding suit- 
able controls. 

Numerical simulations extending the technique to sys- 
tems with non-tridiagonal Hamiltonians suggest that 
constructive control is generally still possible, and sim- 
ilar strong controllability results can almost certainly 

be obtained using very similar arguments for these sys- 
tems. This class would include for example interesting 



systems such as spin-chains with non-nearest neighbour 

couplings. Beyond the extension of generic controllabil- 
ity results to other classes of Hamiltonians. interesting 
questions for future work — the answers to which could 
point the way to achieving effective control with much 
simpler control system designs — include what type of sys- 
tems can be effectively controlled with a single local ac- 
tuator, whether the placement of the actuator matters, 
and how different forms of coupling between the system 
and the actuator affect the control outcomes in practice. 
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